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Abstract. We study the transition to the continuum of an initially bound 
quantum particle in R'*, d = 1,2,3, subjected, for t > 0, to a time periodic 
forcing of arbitrary magnitude. The analysis is carried out for compactly sup- 
ported potentials, satisfying certain auxiliary conditions. It provides complete 
analytic information on the time Laplace transform of the wave function. From 
this, comprehensive time asymptotic properties (Borel summable transseries) 
follow. 

We obtain in particular a criterion for whether the wave function gets fully 
delocalized (complete ionization). This criterion shows that complete ioniza- 
tion is generic and provides a convenient test for particular cases. When sat- 
isfied it implies absence of discrete spectrum and resonances of the associated 
Floquet operator. As an illustration we show that the parametric harmonic 
perturbation of a potential chosen to be any nonzero multiple of the charac- 
teristic function of a measurable compact set has this property. 

Keywords: Ionization, delocalization, resonances, Floquet theory, Borel summa- 
bility. 

1. Introduction 

We consider the non-relativistic Schrodinger equation for the wave fun 

^[-A + V{x)+n{x,t))i, 
where n{x,t) is a time-periodic external potential (not necessarily small): 

n{x, t) = n{x, t + 27r/w), w > 
We take V and Q. real-valued and so that 

VeL°^\M.'^), O e L°°(M'' X [0,27r/w]) 

with ^ satisfying 

f](x, i) = ^ %(a;)e*^'^*, ^■(x) = ViZ~(x)\ sup \^j{x)\f < oo 

We set, without loss of generality, 

D.o{x) = 



^Also Department of Physics. 
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We are interested in the behavior of solutions tp{x,t) for large t when 

^(x,0)=: Vo(2:)€i'(IR''), / |^o|'da; = l (1.6) 

and V'o is sufficiently regular (we assume it of class C^). Of particular interest is 
the survival probability for the particle in a ball B in Mf^, J^^g \tp\'^dx := Vsit)- If 
pB{t) approaches zero as t ^ oo for all B, then we say that the particle escapes to 
infinity and complete ionization occurs. 

While many results in the paper only require (|1.4(l plus sufficient algebraic 
decay of Q and V for large some specific results later in the paper, particularly 
detailed analytic information, require that V and fl are compactly supported, 

supp(y) U supp(f7(-, i)) c D (1.7a) 
and others require the same also for tpo in H1.6|l 

supp(?/'o) U supp(y) U supp(f7(-, i)) C D (1.7b) 

with 

D C M'' compact, M'' \ £» connected, meas(ai:') = 
The rest of the paper will therefore be written in the context of this setting. 

1.1. Nature of the results. Under the assumptions ljl.7b|l '0(a;,i) is obtained for 
large t as a convergent combination of exponentials and Borel summable power 
series in t~^/^. 

If an additional assumption (connected to the absence of discrete spectrum of the 
Floquet operator) is satisfied, the long time expression of ^ contains only decaying 
terms, cf. Theorem 1121 in S 13.81 i.e. we get complete ionization.^ 

We find in Proposition 1141 a convenient sufficient condition for complete ioniza- 
tion, and show that it is satisfied by a nonpcrturbative example^ 

V{x) ^VdXd{x); n{x,t) = 2nDXDix)smujt (1.8) 

where Xd is the characteristic function of D in d = 1,2,3 and Vd and Qjj are 
arbitrary nonzero constants. 

We previously obtained similar results for more general potentials in d = 1 and 
radially symmetric ones in d = 2,3. See [3 |S1 El El E] and |B] where there is a 
review of our previous work on this problem. 

1.2. Strategy of the approach. The key steps of our approach are outlined in 
il3.4l The method we use is based on a study of the analytic properties of V', the 
Laplace transform of ip. The type and position of the singularities of ip^ given in 
Theorem Ul and Lemma jSl provide information about the time behavior of ip; the 
former are obtained from an appropriate equation to which the Fredholm alternative 
approach applies. 



While explicit information on long time behavior requires ■>po to be localized, decay for more 
general ^po G is then an immediate consequence of the unitarity of Schrodinger evolution. 

^The same results hold if a bounded time-independent potential, not necessarily constant, with 
compact support disjoint from D is added to V, see Remark l25l 
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2. Laplace transform, link with Floquet theory 

Exisence of a strongly difFerentiable unitary propagator for Hl.l|l (see [2^1 v. 2, 
Theorem X.71) implies that for tpo e L^(M'*), the Laplace transform 



V'(-,p) / ^P{;t)e~P'dt (2.1) 
Jo 

exists for ^{p) > 0. It satisfies the equation 

(-A + V{x) - ip)^{x,p) = -zV'o - nj{x)tP{x,p~ ijuj) (2.2) 

and the map p —> 4'{'iP) is valued analytic in the right half plane 

p e H {z : 5R(z) > 0} (2.3) 

Clearly, equation H2.2I) couples il^{x,pi) with ijj{x,p2) iff (pi — P2) G iujI^. Setting 

p = i(cr + nw) with 3^(7 e [0,w) (2.4) 

(sometimes it will be technically helpful to relax this restriction on a) we define 
y'n\x; a) = 'ip{x, i{a + nuj)). Eq (|2.2(l now becomes a differential-difference system 

(-A + F + a + nt^)yW = -zV;o - ^ %(a;) (^"'y™)^ (2-5) 

where the shift operator S is given by 

{Sy)n = J/n+l (2.6) 

2.1. Connection with Floquet theory. The solution of Hl.l(l with time periodic 
51 is of course the subject of Floquet theory (see |2] and Cni--23_) and therefore our 
analysis connects to it in a number of ways. Let K be the quasi-energy operator in 
Floquet theory 

{Ku){x, 6) = (^-z^ - A + V{x) + fl(x, e)^ u{x, 6); xeR^,9e (^.7) 
Then, letting 

7/(x,0;a) = ^2/W(x;a)e™"« (2.8) 

riGZ 

be the solution of the eigenvalue equation 

Ku = -au (2.9) 

we get an equation for the ^ which is identical to the homogeneous part of equation 
l|2.5|l'^. Solutions of (|2.9|l with u € L'^(R'^ x S^^^^) correspond to eigenfunctions of 
K. 

Remark 1. If u is an eigenfunction of K corresponding to the eigenvalue —a, then 
^g-'ji^e is an eigenfunction with eigenvalue —a + juj. For this reason it is enough 
to restrict a to the strip given in Jj^.^l ). 



^The functional spaces are different. Proposition^^ clarifies this question. 
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Complete ionization clearly requires the absence of a discrete spectrum of ()2.9|l 
(Otherwise, if u{x, 9) is an eigenfunction of K, then e^'^*-u{x, t) would be a space 
localized solution of the Schrodinger equation.) In a recent work '16', Galtbayar, 
Jensen and Yajima proved that the opposite is also true. They obtained asymptotic 
series in t^^/"^ for the projection of the wave function '4>{x^ t) on the space orthogonal 
to the discrete spectrum of K. 

Our approach via Laplace transform is different from that of |16| . For the cor- 
responding time evolution our results are stronger than those obtained in |16| but 
apply to the more restrictive classes of V and satisfying ljl.7b|) in d = 1, 2, 3. We 
show that the time behavior of il){x^t) is given by a Borel summable transseries 
containing both power law decay and exponential terms. For potentials satisfying 
our condition (|1.8(l we show that K has no discrete spectrum or resonances and all 
the exponentials are decaying. 

More details on the connection between our work, Floquet theory and | 16| are 
given in |S| 

3. Integral equation, compactness and analyticity 

3.1. Laplace space equation. To simplify contour deformation, we first improve 
the decay of ip for large p, by pulling the first two terms in the asymptotic behavior 
for large p from ip. Let <5y = 1 if i = j, and otherwise, 51 j = 1 — 5ij and define 
the operator 91 by 

mUx) = (-A + _ j„,/„_^ + V n^-k{x)^-^^ (3.1) 

a + nio ^—^ a + kuj 

3.2. Green function representation. To pass to an integral form of the system 
of equations (|3.12|) we apply to them the Green function of (—A + a + nuj), given 

by 

(fln/) (x) = J G(At„(x - x'))f{x')dx' (3.2) 

with 
and 



-ip — + mo (when p e H, k„ is in the fourth quadrant) (3.3) 
1 



G(k„x) = < 



^Ko{Kn\x\) d^2 (3.4) 



-^|a;|-ie-''"l"l d = 3 
k 47r' ' 

(see [23) where Kq is the modified Bessel function of second kind. 







Ko{x)^ I e""'°'^'dt^e-=' I ds (3.5) 

Jo Vs(s + 2) 



Note that, in the setting H1.7a|l . for / supported in D we have 

{Snf){x)^ j G{K,,{x-x'))f{x')dx' (3.6) 
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Eq. in integral form becomes 

2/W =-z0„^o + (G:yW)n 

where 



2/™ = (yi^'Wz and (G:2/W)„ 



To ensure better decay with respect to n we further substitute in (|3.8|) 
where (V'i,«)nez V'l and 

V'l = C [(0n^/'o)„ez] 

Then y satisfies 

y = w + £(cr)y 



(3.7) 



(3.8) 

(3.9) 

(3.10) 
(3.11) 

(We write £ for C(cr) when the dependence on a need not be stressed.) In differential 
form H3.1HI reads 



(-A + cr + nuj)yn = iip2,n - Vyn - ^ f^j(a^) (•S ^y)^ 



(3.12) 



3.3. The Hilbert space. To analyze the properties of H3.11|l we use the Hilbert 
space 

n^f^{L\B)) (3.13) 

where B is an arbitrary ball (containing D) defined as the space of sequences 
{Unjnez, yn G i^(-B) with 

and adequate 7; we take for definiteness 7 = 3/2; larger 7 can be taken if one 
assumes more differentiability than (|1.4I) implies. (Note that H is different from 
the Hilbert space L^{L^(R^)) used in Floquet theory. See also fjSl) 

3.4. Strategy of the approach, continued. As mentioned, unitarity of the evo- 
lution shows that tp{-,p) E L'^(R'^) if 3fi(p) > 0. In the integral form (|3.11|) . whose 
solutions are in Ti. when 5R(p) > 0, the operator is, under our assumptions, com- 
pact. The solution of this equation is shown to be unique in Ti. for large enough 
3fi(p) by the contractivity of the integral operator. Uniqueness and analyticity of 
the solution for p in the right half plane H follow by an application of the analytic 
version of the Frcdholm alternative . We then show that the solution and thus 
"0 are analytic with respect to a uniformizing variable, in appropriately chosen do- 
mains containing parts of the imaginary axis. The contour of the inverse Laplace 
transform, t c + ir; r S M, c > 0, can then be deformed to iR (the boundary of 
H) where -0 is analytic except for a discrete set of square root branch points. The 
large time behavior of ip follows. 
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3.5. Compactness. In d 1, wc further transform the equation, see (|8.3|l m 
Appendix to improve the regularity of the operator at n = and a — 0. 

Lemma 2. Under the assumptions <1.7a|) . w €z Ti. and € is a compact operator on 
H. 

To show that w G H we use the fact that the operators g„ satisfy (see Appendix 
A of [Tj, and also JTH . 

sup(l + |n|)i/2||0j|^,(^) <oo (3.14) 

Then 

sup(l + InD^/^IKfln-^o),!!! < oo implying sup(l + |n|)||'i/;i,n||L2(X)) < oo (3.15) 
In view of H1.4|l we also have 



sup(l + \n\) 
nez 



l/Vi,« + ^f^,(a:) {S 



-'y) 



iez 



< oo (3.16) 



implying 

sup(l + |n|3/2)||u;„||^2(p) < oo (3.17) 
nez, 

It is not difficult to check that g„ is compact on L^{B) for each n; it is more 
delicate to show compactness of £; both properties are proved in 13 

3.6. Uniqueness. 

Lemma 3. For large enough — 3(t, eq. I^S.ll]) has a unique solution in Ti. 
The proof is given in Appendix IbI 

3.7. Analytic structure of i/). 

Remark 4. It is convenient to introduce the uniformizing variable a = ; with 
the natural branch of the square root, u is in the fourth quadrant when a is in the 
lower half plane. In this variable, we write kq = u and k„ = V nuj + for n ^ 0. 

Proposition 5. In the setting J1.7a|) the operator £(cr) is analytic in u in the region 
Suj — {u : \'?Rv?'\ < Lu} hence in ^/a in the strip (see Remark^jj 

|cr : 5R((t) e (-w,w)| (3.18) 

Additionally, £(cr) is analytic in a at any a ^ 0. 



Proof. In terms of u we define kq = u, k„ = y'v^+nuj and then k„ is analytic in 
u in the simply connected region Su for any n e Z. Since in the setting l|1.7all the 
integral (|3.6I) is over a compact set, D, g„ is also analytic for any n € Z. Analyticity 
of £ follows from the fact that 



with 



£= lim Cat (3.19) 

and convergence is uniform in u on compact subsets of 5*1^. This is shown in 
Lemma 123 □ 



IONIZATION IN TIME-PERIODIC POTENTIALS 



7 



Proposition 6. There exists a unique solution y to and it has the same 

analyticity properties as € if 

(A) For 5ct < (v ^ €v , v e H) ^ v = (3.21) 

Proof. This is notliing more than the analytic Fredhohn alternative (see e.g. 
Vol 1, Theorem VI.14, pp. 201). □ 

This formulation is convenient in determining the analytic properties of y with 
respect to a, instrumental for the Borel summability results stated in Proposi- 
tion UHl 

Theorem 7. // (A) and fl.7b|) hold, then: 

(i) the solution y of J^8.11]) is meromorphic in u in the disk {u : \u\ < ^/Lu}, see 
Remark^ analytic at u = and in the fourth quadrant of S^. Furthermore, y is 
analytic in a at any (Tq 7^ 0. 

(ii) ^ is analytic in p in a cut neighborhood of iW, {p : ^(p) > —e} with cuts 
toward —00 at incu, n S Z. Furthermore, in a neighborhood of inuj, see 

we have tp{p) = ^n(p) + Bn(j))\/a where An and Bn are analytic at inuj and, for 
some e < ^/uj and \u\ < e we have 

sup n^(\\An{p)\\mB) + \\Bn{p)\\mB)) <<x, (3.22) 

n&.\u\<t ^ ' 

Proof, (i) follows from Propositions [3 and El and from the link between -0 and yn- 
(ii) The functions An and Bn, are simply the even and odd part respectively of 

the analytic function yn{x;u). The estimate follows from the fact that y{x;u) ± 

y{x; —u) G 7i is analytic in u for |m| < e. 

In Proposition 1101 and Theorem 1121 below we use the following result. 

Lemma 8. In the setting <1.7a|) . if ip has a pole at a — ^ iM, then, in the 
variable a — ao if ao =/= or u if ao = 0, the pole is simple. 

This is shown in Appendix lUl □ 

Proposition 9. Condition (A) is satisfied for the potentials 1^1. 
Proof. This is established in f|Hl □ 

3.8. Asymptotic expansion of tp and Borel summability. 

Proposition 10. In the setting fl.7b|) there exist N e N, {Tk}k<N, and {F^;k{t, x)} k 
2-1: / uj -periodic functions of t, such that, for t > Q, 

N 

V'(i,:E) (3.23) 
jel k=i 

with > for all k < N, Pk{t) are polynomials in t, reducing to constants if 
JftFfe = 0, and the hj(t,x) have Borel summable power series in t~^l'^ 

hj{t,x) ^CBY, hkj{x)t-^'^ (3.24) 

with /cq > 1. 
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Remark 11 (Borel summation). If f is a formal power series, say in inverse 
powers of t, then F = Bf is also a formal power series, defined as the term-wise 
inverse Laplace transform in t of f . If (1) F is convergent (2) its sum F can he 
analytically continued along and (3) 3u s.t. F(p) G L^{W'^ ,e~'^Pdp), then the 
Laplace transform CF is by definition the Borel sum of f denoted by CBf. In our 
context we have, more precisely, 

hj {t, x)= Fj (VP, x) e~P* dp hkj {x) t^'^'/^ , t ^ +oo 

•^0 k 

where the functions Fj{s,x) are analytic in s in a neighborhood o/R"*" and for any 
6 €E R there exist a constant C such that for all j and p £ R"^ , 

sup \F,i^,x)e''^\P\\<f, 

p>0;\x\<b 

where the fj decay in j faster than j^^ under the assumption ^1-4}) and factorially 
if Q. is a trigonometric polynomial. Thus the function series in \S.2^) converges 
(rapidly in the latter case). 

The role of condition (A) is described in the foUowing result. 

Theorem 12. (i) If (A) holds, then on the right side of \3.24^ and f.V. g'j)) we have 

fco > 3 and m^u > for all k. (3.25) 

In particular we have complete ionization of the system. (See also Proposition\^ 
as well as Proposition^^ and Remark \15\ below. ) 

(ii) If {A.) is not satisfied, then some ^RFfc may vanish; the part of ^ correspond- 
ing to these Tk remains a spatially localized quasiperiodic function oft. Exception- 
ally, ko = 1 if (A) does not hold (see also |16| and Proposition^^. 

The proofs of Proposition 1 1 01 and Theorem 1121 are sketched in Appendix IdI In one 
dimension a similar result is stated in 0. 



4. Ionization condition for compactly supported potentials 

For the setting H1.7a|l we derive a technically convenient condition implying (A). 
Assume 7^ f G 7i and v = €v. Then there exists a nontrivial solution in Ti. to 
the system 

(-A + cr + nuj)yn = -Vyn - ^ ^j{x)yn-j (4.1) 

We multiply (|4.1(l by integrate over a ball B containing D, sum over n (which 
is legitimate since y G Ti.) and take the imaginary part of the resulting expression. 
Noting that 

j^n^'L j.n^'L j.n^Z, 

= Y ^]{x)y;riym-j (4.2) 
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SO the sum ()4.2|l is real, we get 



n6Z 



ne2 



We take d — Z (the analysis is simpler in one or two dimensions). It is convenient 
to decompose 2/„ using spherical harmonics; we write 

yn- Rn,l,n{r)Yr{eA)- (4-4) 

i>0,|m|<; 

The last integral in H4.3|l . including the prefactor, then equals 

n,m, iRn .771,1 

Rn,m,l] (4.5) 

where rs is the radius of B and VF[/, g] is the Wronskian of / and g. On the other 
hand, since V and are compactly supported, we have outside of B 

Ay„ - (fj + nuj)yn = (4.6) 

and then by 14. 4|) . Rn,i,m satisfy for r > the equation 

, 2 ;(? + 1) 



R" + -i?' 

r 



-i? (cr + ntj)i? 



(4.7) 



where w6 licLVC suppr6SS6d. th.6 subscripts. Let Qn,i,m — rRn,i,m- Then for the gn^i 
we get 

'1(1 + 1) 
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thus 



and 



RR! 



+ (cr + nuj) 



99' \9\ 



5 = 



(4.8) 



(4.9) 



'W\R,R] = W[g,g] =:Wn. 



(4.10) 

Multiplying (|4.8|) by the conjugate of H4.8|l by g and subtracting, we get for 
r > TB, 

W:, = ia-a)\g\^ ^2i\g\^^a (4.11) 

Remark 13. Simple estimates using equation \3.21]) . the definition yj.'^) and 
imply that, for some c„, 

yn{x) = ^-^(cn{e,(t))+0{\x\-^)^ as|x|-^oo (4.12) 



Let us consider two cases of (14.111 . 

Case (i): 3ct < 0. By Remark [T51 we have 



C'e-«"'-(i + o(l)) as r->(X) 



(4.13) 
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There is a one-parameter family of solutions of (|4.8|l satisfying (|4.13|l and the as- 
ymptotic expansion can be differentiated j291 . We assume, to get a contradiction, 
that there exist n for which gn ^ 0- For these n we have, using (|4.13|) . differentia- 
bility of this asymptotic expansion and H3.3|l that 

^ lim \gn\-^Wn = -3k„ > 0. (4.14) 

It follows from (|4.11() and (|4.14() that ^^W^n is strictly positive for all r > and 
all n for which g„ ^ 0. This implies that the last term in (|4.3|l is a sum of positive 
terms which shows that (|4.3ll cannot be satisfied. 

Case (ii): 3(T = 0. For n > there exists only one solution g of H4.8|l which 
decays at infinity (cf. Remark [TSl and the discussion in Case (i)), and since (|4.8() 
has real coefficients this g must be a (constant multiple of a) real function as well; 
therefore we have Wn = for n > 0. 

For n < 0, we use Remark II 31 (and differentiability of the asymptotic expansion 
as in Case (i)) to calculate the Wronskian Wn of 5, g in the limit r ^ 00: Wn — 
|c„p(l -I- 0(1)). Since for Scr = 0, Wn is constant, cf. (|4.11l) . it follows that Wn is 
exactly equal to |c„p. Thus, using H4.3(l and l|4.5|l we have 

yn{x) = for all n < and |a;| > (4-15) 

Proposition 14. In the setting |1.7al) . if (A) fails, then we have 

yn{x) = for all n <0 and x ^ D (4.16) 

Outside D we have Dj/„ — 0, where D is the elliptic operator — A + cr + nu). The 
proof follows immediately from (|4.15|l . by standard unique continuation results |17j . 
|22| . ps ( in fact, D is analytic hypo-elliptic). 

Remark 15. Proposition^^ points toward generic ionization under time periodic 
forcing. Indeed, we see from i4.1()^ that equations i4.1{ l are formally overdetermined 
when n < (i/n is in the domain of A so that, in the function and "one 

derivative" are given on the boundary) and are expected, generically, not to have 
nontrivial solutions even if yn had to satisfy J^S.l^) for n < only. 

The latter reduced problem is relatively easier to study and we used it to show 
that (A) holds in a number of settings, including the potential 1) j.<j^) . see [S]. 

Remark 16. There do in fact exist nongeneric potentials (though not in the class 
Hl.t^) ) for which ionization fails [S] I21L I24j . 

5. Connection with Floquet theory, continued 

Proposition 17. If u is an eigenfunction or resonance^ of the operator K defined 
in {2. 7| ) such that u G L'^{M? x Sl^^i^), then u — €.u inJi and so (A) fails. 

The proof is an immediate consequence of Proposition 1281 in Appendix IeI Con- 
versely, we have the following result. 

Proposition 18. We assume the setting ill.7a|) . 

(i) If V — itv for some do G (0,w) and u G 7i, then v G P(L^(R'^)) thus it is an 
eigenfunction of K . 

(ii) If V ~ £u for ctq = and w G 7i, then, for d = Z, v is of the form v = 
C\x\~^Sno + vix) where i) G P(i^(R'^)) (resonance of K). 



^See 1161 and Proposition 1181 (ii) . 
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Proof of Provosition \lfA For the same reasons as before, we focus on d — 3. 

(i) We see from ():->.6|l . H3.4|l and (|4.16|l that for aU n we have y„ € L'^{ 
Furthermore, a straightforward calculation shows that ||2/n||L2(R3) < C||t/„||i5 where 
C is independent of n. Proposition l28l in Appendix IeI gives the necessary estimates 
in n to complete the proof in this case. 

(ii) For n 7^ we have, for the same reasons as in (i), j/„ G L^(R'^). But now, 
at n = 0, since (Tq = the Green function (|3.4(l does not have enough decay to 
ensure yo g L'^{M.^^). We have instead, for x' E D and |a;| oo, Go{x — x') — 

|a:|~^+0(|a;|"2). The statement now follows from and □ 

6. Example (fO^ 

To show that it can be effectively checked whether (|4.16|l can be nontrivially 
satisfied, we consider the example (ll.8|) . It is convenient to Fourier transform the 
system (14. 1|) in x. In view of (|4.16|1 . for n < 0, y„ = outside D. We then have, 
for n < 0, 



1 

47r 



■„ := / y„e-*'=-^dx = / yne-'^'^'dx (6.1) 



Vn 

and 

-k^Vn = ~e f Vne-^'^-^'dx = / Ay„e-*'=-^rfa: = / Ay^e-'^^'^dx (6.2) 

JR3 JK3 7d 

For the setting (|1.8|l and n < — 1, 14.1|l reads 

{k'^ + (7 + noj)yn = -Voyn + i^D [Vn+l - Vn-l) (6.3) 

Remark 19. For n < —I, the functions y^ are entire of exponential order one; 
more precisely, if B is a ball containing D we have 



\yn{k)\ < VWd) e^"^-- WyJmD) (6.4) 

Proof. This follows immediately from the definition of y. (See also j32j for a 
comprehensive characterization of the Fourier transform of a compactly supported 
distribution.) 

Proposition 20. The generating function 

r(fc,z)= ^i/_„_2(fc)z'" (6.5) 

m>0 

is entire in k and analytic in z for \z\ < 1. 
Proof. Since y E H we have 

||2/«IU^(D) < const |n|-3/2 (6.6) 

Using Remark El the conclusion follows. 

A straightforward calculation shows that Y satisfies the equation 

MY ~ z^ ~ if] ( z - -] Y ^ if3y^i + if3^ (6.7) 
oz \ z J z 

where 

M = Lj~\k^ +a~2Lj + VD) (6.8) 
and f3 — flo/oj. The solution of H6.7|l is 



Y = z'"e 



M-il3(z+z-^) 



(6.9) 
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where the integral fohows a path in which is approached along the negative 
imaginary line. 

Remark 21. Provosition \2(A imvlies C{k) = 0. 

Proof. It is easy to check that otherwise the limit of Y{k, z) as ^ along lEr 
would not exist. □ 
Thus 

y(fc, z) = -z/3z*^e-^(^+-^) £ + ds (6.10) 

We now use the nontrivial monodromy of Y on the Riemann surface of logz, fol- 
lowing from the integral representation H6.10|l . Analytic continuation around the 
origin gives 

)iY{-,ze^-^) y(, z)) = y_i j ds + f ds 

(6.11) 

where C is the curve shown in Fig. 1. Let 

F{M) - i j^ds (6.12) 

Jc * 

Proposition 22. We have 

y-i{k)F{M + 1) + y^2{k)F{M + 2) = (6.13) 
Proof . This follows immediately from the discussion above. □ 

Proposition 23. For every large iV G N, F{z) has exactly one zero of the form 
ZN = N + o(iV°). For large N we have F{1 + zn) ^ 0. 

Proof. It turns out that F{M) is a Bessel function of order M evaluated at 2 
and a proof can be given based on this representation. However, in view of later 
generalizations we prefer to give a more general argument that does not rely on 
explicit representations. 

Let M = N + C, with N a large positive integer, C, complex with |C| = e and e 
positive and small. Let Ci be the counterclockwise circle {z : \z\ — 1} and L the 
segment [0, — i]; we write 

niV + C) = (l-e--<) / — / --j^ds 

Jl 5 Jq^ s 

^27:iC(l + 0{C')) / -J ds+ / -J ds (6.14) 

^ ^ JL S J(j^ S 

where in the integral along L the principal branch of the log is used. The integral 
over L can be estimated with Watson's Lemma, see e.g. 

' ds = jC(27r)3/2/3l-^7V^-3/2g-^g~Cln(/3/A') ^ ^(^0)^ 

By the Riemann-Lebesgue lemma, J^^ ^ as iV — > oo. We get 

F{N + = jA^((2^)3/2/3l-^iV^-3/2g-iVg-Cln(;9/Ar) (^l + o(iV", («)) +o(iV") (6.15) 
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Figure 1. The curve C. 



The existence of a unique simple zero at some N + (^pf with jC^v I < e is a consequence 
of the argument principle. 

The position of ^jv can be found more accurately as follows. We have 



c 



N 



-ds 



-ds] 



(6.16) 



from which it follows that (n ~ o{ const^/iV!) which readily implies that (|6.16|l is 
contractive and that 



-ds 



pi/3(s+s 



-ds 



1 



o(iV°)) 



(6.17) 



Using (|6.15|l and the fact that the first integral in H6.17|l gives the Laurent coeffi- 
cients of 6*'^'-''+'* ) which can be independently estimated from the series expansion, 
we find that, with constants that can be calculated. 



|C^|=cic^iV-2^+«(l + o(7V«)) 



(6.18) 



Thus (n+i/(n — > as —> oo and the second part of the Proposition follows. 
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Proposition 24. Relation Iff, j'jl) . with y-i{k),y-2{k) entire of exponential order 
one (cf. Remark \lfJ\) implies 

y-i(fc) =t/-2(fc) = Vfce (6.19) 

and then 

yn{x) = Vn e Z and almost all x e M.'' (6.20) 

Proof. Propositions 1221 and 1231 and H6.8|l imply that y^i{k) has at least const 
zeros in a disk of large radius R. Since y_i(fc) is an entire function of exponential 
order one, it follows that y^i = (see, e.g. [3S]). By (|6.13|) we have y_2 = 0, so 
that y_i = y_2 = . 

In the present model (|4.1() reads 

(-A + o- + nuj)yn = -VoXoyn - i^DXaiyn+i - Un-i) (6.21) 

and H6.21|l with n = —1 and n = —2 implies XdUq = 0; XdVs = respectively; 
inductively XDUn — for all n. Then (—A + a + nu!)yn = in R'' \ dD. Since, by 
Proposition UHl y„ G L'^{R^) we see (for instance by taking the Fourier transform) 
that ||y„|| ^ 0. 

Remark 25. R can be shown that condition (A) holds with VdXd replaced by 
VdXd +^i(a^) where Vi{x) is bounded, not necessarily constant, with compact sup- 
port disjoint from D. On the support of Vi, is zero and it can be seen that 
for n sufficiently negative yn is zero on the support of Vi . From this point on the 
arguments are very similar, but we will not pursue this here. 



7. Compactness 

7.1. Compactness of the operator g„ defined in (|3.6|) . The case d = 1 is 
discussed in Appendix 1X1 For c? = 2,3 compactness follows from Theorem VI. 23 , 
Vol. 1, pp. 210 of |S| (for d = 3, note that e-'^^^'^-y^ /\x - y\ & L^{D x D)). 

7.2. Compactness of £. The property (|3.14(l is mentioned in Appendix A of ^J. 
We include here an elementary proof of (|7.1(l below (which also can be refined 
without serious difficulty to yield the sharper estimate (|3.14f) '). 

Lemma 26. We have 

llflnll — >0 as \n\ ^ oo (7.1) 

(where \\ ■ \\ is the L^{D) ^ L^(D) operator norm) uniformly in u in the region 
Suj,e,A — {u ■ \u\ < A, |5R(m^)| < ll) — e}, where A > and e is any small positive 
number. 

Proof. Relation (|7.1|l follows from a general result by Agmon 01 which provides 
estimates on the rate of convergence. We give below an elementary proof in our 
case. 

We prove the result for d — 3 (the proof is simpler in d = 1, 2, noting that for 
large x with arg x € (— 7r,7r) we have Ko{x) = yj^e~^x^^^^{l + o(1))). Define 

Q^{x',x")^ dx—-^ — (7.2) 
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We have 



||g„f = sup / Qn{x\x")f{x')fix")dx'dx" 

11/11=1 •''c^ 

< (^J JQnix\x")\^dx'dx'^ (7.3) 
The last integral goes to zero as |n| oo. To see that, note that 



y/ nuj + u2 = +0{n ^/^); as n ^ +oo (7.4) 
and using the triangle inequality we get 

dx 



|Qn(x',x")| < Const e~^l^'-^"l 



\x' — x\ \x 



< Const 6-^1^'-=""! (7.5) 



and the conclusion, for n — > +00 follows by dominated convergence. We now focus 
on large negative n. Since 



■\/ nuj + = —iy^\n\uj + 0{n ^^^); as 71 — )■ — 00 (7.6) 

it is easy to check that (|7.1|) follows once we show that ||£|[^]|j —> as cxd where 
is obtained by replacing k„ with iv in the definition of g„. We first show, 
with an analogous definition of that 

sup \Q[^]{x,x')\ = Qa < 00 (7.7) 

Indeed, we choose a ball Bi, centered at x' of radius b large enough so that it contains 
D and write the integrals H7.2() in spherical coordinates centered at x' with x" on 
the z axis; in these coordinates |a; ~ x'\ — r and \x — x"\ > d{x, Oz) — rsmO and 
thus 

\Q[.]{x,x')\< f ^ -<f drd0dcf> < Anb (7.8) 

JBt \X - X \ \X - X I Jb^ 

Let p{x; x' ,x") ^ \x ~ x'\ — \x — x"\. We then have \p{x; x' ,x")\ < \x' — x"\ and 
we get 

r pii'p{x]x' ,x") r\x ' —x \ 

Qy,^(x\x")= \ dx- — 1= / e-''dM(p) 

Jd \x'-x\\x"-x\ J_\^n_^,\ 

where the positive measure [i is defined by 

/i(A) = M.',."(A) = / f^^^ : (7.9) 

J {x:p(x)eA)ir\D \x — x\\x — x\ 

Since the integrand in (|7.9|1 is in L^, the measure [i is absolutely continuous with 
respect to the Lebesgue measure m. We let /i(p; x', x") = then h ^ and we 
get 

Q[,](.t',x")= /' ^ e"'Ph{p;x\x")dp (7.10) 

J-\x"-x'\ 
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By the Riemann-Lebesgue lemma we have^ 

Q\^y]^{x' ,x") ^ as — > oo (7-11) 

Now H7.7(l . (|7.1H) and again dominated convergence implies \\Q[u]\\ — > as u — > oo 
completing the proof. 

Lemma 27. Under the assumption fl.7a|) . the operator £ is compact on Ti. and 
analytic in u in S^^^^.A, cf. Lemma W^ 

Indeed, € is the norm limit (|3.19|l . uniform in u e Suj,t,A, where <Ln are compact 
by Lemma ESI and analytic as explained in the proof of Lemma We note that 
the operator 



is bounded in Ti.. Indeed, if we write (n) := 1 
{n) < (j) {n - j) and 

2 



(7.12) 

|n| we have, for (n,j) G 



E 



< 



1 



(ji> ^ I O2) 2 1% I E " -^'i) ' I " -^s) ' l2/"-J2 1 



<bllp|IOT/'%ll?<C'||y||p (7.13) 



by p.4|l . □ 



8. Appendixes 



A. Compact operator formulation and proof of Lemma for d = 1. We 

can assume without loss of generality D C [—1,1]. For n = we choose some large 
positive a such that sin2Y^ 7^ 0, denote by f±{x) the functions e^"^ and let tp^- 
be the solution of the equation 

-7/;" + (aX[_i4]+u2)V = (8.1) 

(see Remark 0J with initial condition V'+(l) = /+(l)i '0+(l) = f'+i^)^ and similarly 
let ip- be the solution of (|8.1|l with initial condition ?/;_(— 1) = /_(— 1), jp'_{—l) = 
/_'_(— 1). Since both the equation and the initial conditions are analytic in u at 
u = 0, so are the solutions 4'± and their Wronskian W{u). It can be checked that 
W{0) — Y^sin2Y^ 7^ 0. In fact, taking Tq = \/ a — we have 

In a neighborhood of u = we write for n 7^ the same integral expression H3.11|l . 
while for n = we write 



(8.3) 



^Noting that the estimate of the norm of 0„ can only increase if extended to (B) where B 
is a ball containing D, and that h calculated in B is piecewise smooth the max of Q[^] can be in 
fact bounded by an inverse power of u; we do not however need this refinement here. 
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where 



(8.4) 



With the same conventions, we now write the integral system in the form H3.11|l . 
Compactness and analyticity are now shown in the same way as for d = 2, 3. □. 



B. Proof of Lemma |3j Let a = (Jq — 2iT where ctq G [0,uj) and r > 0. We show 
that 

||e:|| ^0 as T^oo (8.5) 

and uniqueness foUows by contractivity. The calculation leading to (|8.5|1 is quite 
straightforward, but we provide it for convenience. In d = 1,2 the estimate follows 
from the behavior of the Green function for large argument. We then focus on 
d = 3. By H3.3|l we have 



Q (('^o + nujf + 0-0+ nuj)^^'^ + 



1/2 



(8.6) 



For ri > we then have SR^n > y/nuj and the same calculation as for (|7.5(l shows 
that 

as n — > +00 



7) 



uniformly in r. For n < H8.6(l gives 

\Qn{x',x")\ < \Q.{X',X")\ 

where 

—i> :~ 3(Kri) ^ oo as n — > — oo 
and now (|7.11|) shows that 

llflnll — > as n ^ -co 

uniformly in r. We choose then no large enough so that 

sup ||0„|| < e 

n>no ,T>0 

For T large enough we have, still from (|8.6|) . 

Choosing a ball B centered at x containing D, we then have for large r and some 
constants independent of /, t, x and n G (— no,no), with the notation a = ^t^/^, 



(8.9) 
(8.10) 

(8.11) 
(8.12) 



(fin/) (a 



< 



-a\x—x'\ 



\\f\\L^{D) 

< C'i||/IIl2(_d) 



<^\\f\\LHD)<e (8.13) 



and the conclusion follows. 
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C. Proof of LemmalSl Assume ao is a value of a where invertibility of / — C((To) 
fails. Then do G [0,cj). By the Fredholm alternative we know that / — £(ct) is 
invertible in some punctured neighborhood of (Tq where the solution of (|3.11|) is 
meromorphic. 

(i) : (Tq 7^ 0. Denote ( = a — <jq. We rewrite H3.11|l in a suitable way near ctq. We 
have from H3.12|l 

(-A + CTo + nuj)yn = -i^2,n - Vyn -CVn + Y^ ^j{x) (S'^y)^ (8.14) 
which we write symbolically 

Wy = -Cy - i^P2,n (8.15) 
and from and |TTT|) we have ************* 

Vn = -i3n1p2,n " Qn Vy^ " C^n " ^ {^~^V)n {^-IQ) 

implying the following version of l|3.11|l . with evident notation, 

2/ - yo + Cfly + ^{<jQ)y (8.17) 

On the other hand, 

oo 

y= c,C^" (8.18) 

with the coefficients Cj G H. Assume, to get a contradiction, that M > 2. Inserting 
in (|gT7|) we get 

c_M - C:(fTo)c-M (8.19) 

C-M+l = C(cro)c_M+l + QC-M 

In differential form we have, 

Wc-M = (8.20) 

Wc^M+i = -c-M (8.21) 

By ifSTT^ and Proposition [THl we have c_m G ^^(^^(M'*)) := Hi (in fact, (c_m)« 
decay at least exponentially in |a;|). On the other hand we then have from (|8.2(J|) 
and noting the formal self-adjointness of 22J, 

(c_a/,c_m) = -(^c-M,Wc^M+i) = -(Wc-M,c^M+i) = (8.22) 

which is a contradiction. 

(ii) (To = 0: there are two differences w.r.t case (i): (a) meromorphicity and 
Laurent expansions now use the variable u — ^/ct; and (b) c_m is not necessarily 
in Hi so we work with Hb — 1^{L'^{B)) for large enough B. These differences can 
be dealt with straightforwardly, so we just outline the main steps. The Laplacian 
is the only ingredient of 20 not formally self-adjoint in Tis. Integration by parts, 
implicit in H8.22|l brings in boundary terms of the form 



fVg ■ dS (8.23) 

where / and g are c_m or c_m+i. Both / and g have decay and this behavior 
is differentiable, as is manifest from H8.19|l . 1)3. 4|l . and 1)3. The contribution 
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from the integrals l|8.23|l is thus 0{rg^) which equals the norm ||c_a/||-Hb, clearly 
nondecreasing in rs- This again forces c_m = 0, a contradiction. □ 

D. Sketch of the proof of Proposition [Tol and Theorem I12L We first show 
Theorem 1121 (i). The contour of the inverse Laplace transform can be deformed as 
shown in Fig. 2. Pushing the contour of integration to the left brings in residues 
due to the meromorphic integrand, and since the kernel of the inverse Laplace 
transform is e^', residues in the left half plane give rise to decaying exponentials 
in 'ip{x,t). Uniform bounds on the Green function as p — > —oo are easy to prove. 
Consequently, there are only finitely many arrays of poles of tp. The contour of 
integration in the inverse Laplace transform can be pushed all the way to —oo in 
view of the exponential decay of the kernel e^*. We are left with integrals along 
the sides of the cuts which, after the change of variable p —p (or p <-> — pe*" if 
poles exist on the cuts), are seen to be Laplace transforms. Since ^ is analytic in 
^/p + into, the contour deformation result shows, ipso facto, Borel summability of 
the asymptotic series of •ip{x,t) for large t. 

The general case is proved in a very similar way, using Lemma |H1 If (A) does 
not hold, then some of the poles of the meromorphic function yn{<^) can be on the 
segment (— w, w). If a pole is placed at cr = 0, then analyticity in u in the operator 
entails a singularity of the form <7^^^^A{<7) + B{a) with A and B analytic, whence 
the conclusion. 

E. Estimates needed for Proposition [T7I and Proposition [T51 

Proposition 28. Let y be a solution in P(L'^{B)) of the homogeneous equation 
associated to \3.1'£\] . Under the assumptions fl.7b|) . we have 

bj||L2(i3) =0(r') as\j\^^ (8.24) 

Let e be small enough and choose jo > large enough (the proof is similar for 
jo < 0) so that ||flj||L2(B) < e for all j > jo, see (|7.1|) . We consider the Banach 
space Bjg of sequences of functions {yj}j>j„ defined on B for which the norm 

\\y\\j„ := sup f\\yj\\L2(^B) (8.25) 
j>jo 

is finite. For j > jo we write the homogeneous part of H3.12|l in the form 

m>0 

+ ^iyj-i{x)+ Y ^iVj-ii^) 

= SjVyj + iSj [ Y ^-rnyrn+]ix) + ^ Vtiyj^i{x) + Ej{x) 
m>0 0<l<j-jo 

=:(%y + £;)^ (8.26) 

Since ||rij||i2(3) — 0{j^^) and y € P{L^{B)) we see that \\E\\j„ < oo. It can be 
checked that 1j„ : Bjg Bjg is bounded, that ||Tjo|| — > as jo oo, and thus eq. 
(|8.26|l is contractive if jo is large. The Proposition follows. 
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and suggestions, including the argument in H7.13|l . We thank A. Soffer for helpful 
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Figure 2. Deformation of the integration contour. 
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Welcome to the deleq package! 



This is a short document to demonstrate the use of the deleq package and its commands. 
It uses deleq version 4.41 (July 7, 1997). deleq was written by Mats Dahlgren (matsd@sssk.se 



http://www.homenet . se/matsd/). Suggestions for improvements and bug reports are most 



welcome, see the documentation, deleq is fully compatible with the leqno option and most 
of the f leqn option. 

We start this demonstration by a simple and well-known equation to get the equation 
number counter going: 

sin^ a + cos^ a = 1 (1) 

The first example will be to make use of the deqn environment to get a partially numbered 
equation: 

sin(— a) = — sin a (2a) 
This equation has its cosine companion, here written in the ddeqn environment: 

cos(— q) = cos a (2b) 

Not so exciting, so far. :-) 

In the next example we introduce the deqarr environment for writing equations: 

sin(a + 13) = sin a cos /3 + sin /? cos a (3a) 
sm{a — P) = sin Q cos /? — sin /? cos a (3b) 

where we also have put in the label Demol in the second equation. Next, notice how the 
environment ddeqar uses the same main equation number as the previous equations: 

cos(q + /?) = cos a cos /? — sin Q sin /3 (3c) 
cos(q — /?) = cos a cos /? + sin a sin /3 (3d) 

In the first of these equations, an \arrlabel command with the label Demo2 is included. 
Now, we will make another eqnarray-like structure, again in the ddeqar environment: 

sin 2a = 2 sin a cos a (4) 

cos 2a = cos^ a — sin^ a (5a) 

or 

= 2 cos^ a - 1 (5b) 

This example shows the use of \heqno in the first equation, which produces an ordinary 
equation number. The first equation is also labelled, with the label Demo3. The next equa- 
tion's number was produced by the command \nydeqno. Also, notice the use of \rem{or}, 
which results in the text "or" without changing the alignment. To illustrate the use of 
\arrlabel{Demo2} above, we here make a reference to it: ... in equations 01 • . .which was 
created by typing '\ldots-C]- in equations \ref {Demo2]- \ldots'. 

To step the equation number counter, we want the following equations typeset in ordinary 
eqnarray environment: 



2 sm a 

tan a = ^ — (b) 

cos^ a 

sin a 

tana = (7) 

cos a 



1 



The upper equation was given the label Deiiio4. 

Now one of the "recycling" commands is to be demonstrated. First of all, let's make an 
ordinary reference to equation |3 and then "recycle" it: 

sin 2a = 2 sin a cos a (j^ 

This was obtained with $$ ... $$ with the command \reqno{Demo3} at the end. Also par- 
tially numbered equations can be recycled with the \reqno{FOO} command, as with equa- 
tion |^ 

sin(a — (3) = sin a cos P — sin (3 cos a ()3b|) 

(However, the use of \rndeqno{FOO} and \rdeqno{FOO} will produce strange results with 
two (different) partial equation numbers if FOO refers to a partially numbered equation.) 
The next equation was written with $$ ... $$ and a \deleqno command at the end: 

2 tana 

tan 2a = k — (7a) 

1 — tan^ a 

Notice how the main equation number counter is still the same, despite that we now are 
outside of the ddeqar environment. Also, notice how the use of \reqno above did not affect 
the equation number counter. The following equation is also set within $$ ... $$, but it 
uses the command \nydeleqno to produce an equation number with a new main number: 

sin^ a = 1 — cos^ a (8a) 

Now we will elaborate a little on \riideqno{FOO} and \rdeqno{FOO}. Let us use equa- 
tion IHl and see what happens if we use the \rndeqno{FOO} command: 

2 sin^ a — 

tan a = — (|6fi) 

cos^ a 

which we rewrite: 

2 sin^ a — 

tan a = ^ — ISp) 

1 — sin a 



with \rdeqno{Demo4} at the end. This is great fun, so why not one more: 

.2 

cos^ a 



2 1 — cos a ^ , 

tan a = ^ &> 



Also, the middle form of equation got a label, Demo5. In the last example, the construct 

\ddeqreqno [-\ j otbaseline] {Demo4} \nonumber 
in a deqrarr environment is used, to obtain f leqn compatibility. 

There is still one equation to write. This time we again use the deqarr environment: 

cos a 

cot a = (9a) 

sm a 

where 

= (9b) 
tana 

This example also shows the use of the \where command, which is a special case of the \rem 
command. Now, the interesting thing of referring to the recycled equations is ahead. If one 
writes \ref -[Demo5}, this is what E^T[^ will return: That is not too instructive, since 
the partial equation numbers are quite common in this document. To make the complete 
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reference, use \ref {Demo4}\ref {Demo5}. The reference 00 is much more comprehensible, 
right? (If you get bad hne-breaks at such references, put them in an \mbox{. . .}.) 

Now we will show the use of the commands to recycle equation numbers in eqnarray-like 
structures. These commands are \eqreqno{FOO}, \deqreqno{FOO}, and \ddeqreqno{FOO}. 
We start by repeating equation l,3bl in an deqrarr environment: 

sm{a — P) = sin a cos /? — sin /? cos a ()3b|l 
sin 2a = 2 sin a cos a 

The second equation here shows how the alignment of the equations is preserved. Next, we 
can use the equation 0] for some variations: 

sin 2a = 2 sin a cos a (^t-) 
sin 4a = 2 sin 2a cos 2a 

= 2 ^2 sin a cos a ^cos^ a — sin^ ^p) 

= 4 ^sin a cos a (^cos'^ a — sin^ a^ ^ ^) 

The features used here are \deqreqno{Demo3} on the first line, \ddeqreqno{Demo3} on the 
second line, and \deqreqno [-\jotbaseline] {Demo3} on the third. Note how the inclusion 
of [-\jotbaseline] on the third line prevents I^T^5]X from inserting an extra blank line. The 
second line ends with '\\', and has thus no equation number at all. 

The last part shows how you may use \eqreqno and its relatives in deqarr environment: 



sina = sin Q cos + sin cos a (10a) 

sin 2q = 2 sin a cos a @ 

sin 3a = sin 2a cos a + sin a cos 2a (lOb) 

sin 4a = 2 sin 2a cos 2a (10c) 

sin 5a = sin 3a cos 2a + sin 2a cos 3a (@ll) 

sin 6a = 2 sin 3a cos 3a (lOd) 



The interesting part here is an '\eqreqno{Denio3}' at the end of the second line and a 
'\ddeqreqno{Denio3}' at the end of the fifth line. Notice how the partial equation num- 
ber counter in unaffected by the enetering of the deqarr environment. 

Finally, the £J\fT>\ If you want your equation numbers left-aligned, just specify the leqno 
option for the documentclass you are using. It should work! If you have any suggestions, 
corrections or contributions, please contact me. Enjoy I^T^5]X! 

mats d. 
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